A finite range pairing force for density functional theory in superfluid nuclei 
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The problem of pairing in the 1 So channel of finite nuclei is revisited. In nuclear matter forces of 
separable form can be adjusted to the bare nuclear force, to any phenomenological pairing interaction 
such as the Gogny force or to exact solutions of the gap equation. In finite nuclei, because of 
translational invariance, such forces are no longer separable. Using well known techniques of Talmi 
and Moshinsky we expand the matrix elements in a series of separable terms, which converges quickly 
preserving translational invariance and finite range. In this way the complicated problem of a cut-off 
at large momenta or energies inherent in other separable or zero range pairing forces is avoided. 
Applications in the framework of the relativistic Hartree Bogoliubov approach show that the pairing 
properties are depicted on almost the same footing as by the original pairing interaction not only 
in nuclear matter, but also in finite nuclei. This simple separable force can be easily applied for the 
investigation of pairing properties in nuclei far from stability as well as for further investigations 
going beyond mean field theory 
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Along with improved techniques to investigate more 
precisely nuclear systems being considered as well known, 
the recent generation of radioactive beam facilities en- 
ables us to examine exotic systems with extreme isospin 
values. Therefore experimental and theoretical studies 
of nuclei far from the valley of /3-stability are presently 
at the forefront of nuclear science. Experiments with ra- 
dioactive nuclear beams have already in the past discov- 
ered a number of new structure phenomena in exotic nu- 
clei with extreme isospin values, and the next radioactive- 
beam facilities in construction will present new exciting 
opportunities for the study of the nuclear many-body sys- 
tems. 

Nuclei far from stability have also an important in- 
fluence on astrophysical processes. Therefore the study 
of such nuclei has wide ranging applications in modern 
nuclear astrophysics. Unfortunately many of the nuclei 
of interest in this context have such large neutron ex- 
cess, that it will be impossible in near future and prob- 
ably even excluded in far future to investigate them on 
earth by experiments in the laboratory. Therefore it is 
extremely important to provide a powerful theory for a 
reliable description of nuclei close to the limits of sta- 
bility. It should be based on a consistent treatment of 
both ground and excited states and should allow for pre- 
dictions of nuclear properties in areas, which are hard 
or impossible to access by future experiments. Ab initio 
calculations and multi-configuration mixing within the 
shell- model are definitely a goal, but, so far, they can 
only be applied in light nuclei. At present, for a univer- 
sal description of nuclei all over the periodic table, Den- 
sity Functional Theory (DFT) based on the mean-field 
concept provides a very reasonable concept. DFT has 
been introduced in the sixties in atomic and molecular 
physics 0, Q and shortly after that in nuclear physics 
under the name 'density dependent Hartree-Fock the- 
ory' 0, 0] ■ Today it is widely used for all kinds of quan- 



tum mechanical many-body systems. DFT can, in princi- 
ple, provide an exact description of many-body dynamics, 
if the exact density functional is known, but for systems 
such as nuclei one is far from a microscopic derivation 
and the most successful applications determine the func- 
tional in a phenomenological way. Starting from basic 
symmetries the parameters are adjusted to characteris- 
tic experimental data in finite nuclei and nuclear matter. 
In addition nuclei are self bound systems. One usually 
considers densities in intrinsic frames and it is still under 
debate whether density functional theory can be exact 
under these circumstances 0, @ . Nonetheless in practice 
DFT provides in many nuclei all over the periodic table 
an amazingly successful description of the complicated 
many-body system [7|, [8( . 

Conventional DFT with a functional E[p] depending 
only on the single particle density can be applied in nu- 
clear physics practically only in a few doubly closed shell 
nuclei. In all nuclei with open shells, and this is the vast 
majority of all nuclei in the periodic table, the inclusion 
of particle-particle (pp) correlations is essential for a cor- 
rect description of structure phenomena. Although, in 
principle the effective pp-interaction is isospin dependent 
with aT = and a T = 1 part, for the vast majority 
of pairing effects in nuclei only the T = 1 part is im- 
portant. In fact, little is known on the effective T = 
part. It is still an open question, whether the effective 
pp-interaction in the T — channel is stron g en ough to 
produce a pairing condensate with We 
therefore restrict ourselves in the following discussion to 
T = 1 pairing correlations between like particles. 

In the framework of DFT pairing correlations are ta- 
ken into account in the form of Hartree-Bogoliubov the- 
ory, where the energy functional E[p, n] depends not only 
on the normal density p = (a + a) but also on the pair- 
ing density k = (a + a + ). Both densities can be combined 
to the so-called Valatin density 1Z [l4[ with the property 
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7Z 2 = 7Z. This shows clearly that the Hartree-Bogoliubov 
version of DFT is a generalized mean field theory. In 
Ref. [ll d El relativistic Hartree-Bogoliubov (RHB) 
theory has been introduced for the treatment of pairing 
correlations in relativistic DFT. It turns out that pair- 
ing itself is a non-relativistic effect influencing only the 
vicinity of the Fermi surface. The effect of the pairing 
field on the small components can be neglected to a very 
good approximation [18J. 

For nuclei close to the /^-stabilit y li ne, pairing has been 
included in non- relativistic DFT 191 or in the relativis- 
tic DFT [20( in the form of the simple constant gap ap- 
proximation, where the pairing gap A is obtained from 
odd-even mass differences. The occupation numbers v\ 
are then determined by the BCS-ansatz. The pairing 
tensor is diagonal with the elements — UfeVk where 
u 'k + v k = 1 an d the pairing part of the density func- 
tional is given by the trace of k: E pa i r [K] = — A^ fe 
Of course, this sum diverges and therefore one has to 
restrict the sum over A: to a pairing window. This intro- 
duces an additional parameter, which is not well deter- 
mined by experiment. Several prescriptions can be found 
in the literature to treat the pairing window [2l|, [22] • 

This way to include pairing correlations corresponds to 
the seniority model [23|], which is also sometimes called 
monopole pairing, where only pairs coupled to angular 
momentum J = feel the effective pp-interaction. Of 
course, this is a very simplified description of pairing 
correlations in nuclei and therefore this force is often re- 
placed by a zero range force 

V^ = V S(r x -r 2 )L(l-P-) (1) 

which is sometimes chosen to be density dependent [24j . 
For zero range forces as in Eq. (1) the factor |(1 - P a ) 
projects on the 1 Sq channel. This force is simple to 
handle in r-space and it is therefore often used in non- 
relativistic HF-BCS-calculations with Skyrme forces [25| 
and in relativistic density functionals based on point- 
coupling models [2|| H3, H^]. Unfortunately the corre- 
sponding pairing energy diverges in this case too. As 
in the case of seniority pairing one needs a pairing win- 
dow. Sharp pairing windows have the tendency to lead to 
flip-flop solutions in self-consistent theories and therefore 
one uses in most of the present applications soft pairing 
windows [29l | with rather arbitrary cut off parameters. 

For nuclei far from stability the BCS approximation 
presents only a poor approximation. In particular, in 
drip- line nuclei the Fermi level is found close to the parti- 
cle continuum. The lowest particle-hole (ph) or particle- 
particle (pp) modes are often embedded in the contin- 
uum, and the coupling between bound and continuum 
states has to be taken into account explicitly. In these 
cases the BCS model does not provide a correct descrip- 
tion of the scattering of nucleonic pairs from bound states 
to the positive energy continuum because several levels in 
the continuum become partially occupied leading to a gas 
of nucleons surrounding the nucleus. Including the sys- 
tem in a box of finite size leads to unreliable predictions 



for nuclear radii depending on the size of this box. In the 
non-relativistic case, it has been shown that the Hartree- 
Fock-Bogoliubov (HFB) theory in the continuum pro- 
vides a very elegant solution to this problem [3(1 HH, [Hj] 
and this method has been applied also to investigations 
of drip line and halo nuclei within covariant DFT [33l.[3^|. 
The HFB theor y pr esents a unified description of ph- and 
pp-correlations [351 ] on a mean-field level by using two av- 
erage potentials: the self-consistent Hartree-Fock field T 
which encloses all the long range p/i-correlations, and a 
pairing field A which sums up the pp-correlations. 

In order to avoid the complicated problems of a pair- 
ing cut-off Gogny derived his energy functional [3] from a 
finite range force of Brink-Booker type. The finite range 
guarantees that the force decreases as a function of the 
momentum transfer and the gap equation converges with- 
out any problems. Thus a pairing cut off is not necessary. 
The parameters of this force have been adjusted very 
carefully in a semi-phenomenological way by Gogny [j] 
and his collaborators [H, [37| to characteristic properties 
of the microscopic effective interactions and to experi- 
mental data. Over the years this method has turned out 
to be a very successful way to describe pairing correla- 
tions in nuclei and it is often used as benchmark for more 
microscopic investigations (38t [39j . 

Of course, mean field calculations with finite range 
forces turn out to require a substantial numerical effort, 
in particular in three-dimensional applications to triaxial 
nuclei and to rotating systems [U [43[ or in ap- 
plications going beyond mean field such as the Generator 
Coordinate Method (GCM) [44[ connected with projec- 
tion to goo d angular momentum [27l . |45| and particle 
number |28l . |46| . Since one needs nowadays systematic 
investigations over a wide range of nuclei [47| it is highly 
desirable to find an effective interaction in the pairing 
channel which is numerically simpler without loosing the 
nice properties of this force. In this investigations we 
propose a new realistic pairing force, which is carefully 
adjusted in a separable form in momentum space to nu- 
clear matter properties of the Gogny force. It turns out 
that this force is simple enough, that its matrix elements 
in finite nuclei can be expressed as a rather limited sum 
of separable terms. 

We start our investigations in symmetric nuclear mat- 
ter. The effective pairing force in the pp-channcl can be 
represented by a sum of all diagrams irreducible in pp- 
direction. In lowest order it corresponds to the bare NN 
interaction. Recently Duguet [48] proposed a microscopic 
effective interaction to treat pairing correlations in nu- 
clear matter in the channel. Starting from an ansatz 
separable in momentum space, he derived after several 
approximations an effective pairing force with zero range 
for practical applications in the context of Skyrme cal- 
culations in r-space for finite nuclei. We start from a 
similar ansatz in nuclear matter, however, instead of re- 
ducing it to zero range, we transform the force obtained 
in this way from momentum space to r-space. Then cal- 
culations in finite nuclei are carried out in terms of an 
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expansion in the eigenfunctions of a harmonic oscillator. 
This is particularly simple for a Gaussian ansatz leading 
to a naly tical expressions for the matrix elements. In Refs. 
[49l |50| | a similar technique was used based on spherical 
Bessel functions. 

The wave functions in infinite nuclear matter are char- 
acterized by the momentum k and the spin s 



Vk s (r) 
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The antisymmetric matrix element of the pairing inter- 
action V se p° in the plane-wave basis has the form 



(k lSl k 2 s 2 |i(i - p a )v s °|k' lS ik 2 4) a 
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where 



(3) 
(4) 



(k^lV 13 "^) = (k|l/ s °|k')(2^) 3 <S(K - K'). (5) 

K= ki + k 2 and k =g(k 1 — k 2 ) are the total and relative 
momentum of a particle pair, respectively. It has been 
pointed that the bare interaction in the 1 Sq channel is to 
a good approximation separable and nonlocal at low en- 
ergy (5lj | . The center of mass part of the matrix element 
is therefore approximated by a separable form, 



(ki^/iik') = -Gp(k) P (k') 



(6) 



The isospin quantum number is not specified in this ex- 
pression, since the form of the matrix elements in the 
(T = 1) channel is trivial. G is the strength parameter 
of this pairing interaction. 

In the 1 S'o channel we find the following gap equation 
in the plane wave basis, the usual BCS equation (52l.l53l]. 



A(fc) = - / 
Jo 
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with the quasi-particle energy 



E(k) = V(e(fc)-^) 2 + A 2 (fc) 



(7) 



(8) 



and the in medium on-shcll single particle energies e(fc) 
associated with the state tp^. /i is the chemical potential 
determined by the density and v{k,k') = —Gp(k)p(k'). 
For a given pairing interaction, one can solve the BCS 
gap equation and calculate the corresponding gaps as a 
function of the density, i.e. as a function of the Fermi 
momentum in nuclear matter. This relationship between 
the gap at the Fermi surface and the Fermi momentum 
determines the properties of the pairing correlations. 

Inserting the separable interaction Gp{k)p{k') into 
Eq. ([7]), the solution of the gap equation is trivial A(/c) = 
Aop(fc), where Ao is the gap at zero momentum satisfy- 
ing the equation 
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V(e(fc)-M) 2 + Agp 2 (fc) 



(9) 
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FIG. 1: Comparison of 1 So pairing gaps at the Fermi sur- 
face A(k F ) as a function of the density for the Gogny forces 
Dl and D1S (solid curves) and the corresponding separable 
forces (dashed curves), where the function p(k) in Eq. (|6]) is 
represented in panel (a) by one Gaussian in Eq. (|15[1 and in 
panel 3 by a sum of three Gaussians in Eq. (I13p . 



and therefore depending on the Fermi momentum k F - 
The gap at the Fermi surface is obtained through 



A(k F ) = A (k F )p(k F ) 



(10) 



which shows that the bell shaped curve A(k F ) determines 
the form of the functional dependence of the function 
p(k) in Eq. @. 

For RHB calculations in finite nuclei it is common to 
use the Gogny force or a zero-range ^-interaction 
[26| . We apply this method to derive a separable forces, 
i.e. the function p(k) in Eq. ^ by mimicking the non- 
relativistic Gogny force in the channel of nuclear mat- 
ter. Of course this procedure depends on the self energies 
e(k). When we deal with the Gogny force this self ener- 
gies Eqs. ([7]) and are no longer free single particle 
energies. Medium corrections have to be included, using 
for instance the Brueckner-Hartree-Fock approximation. 
For simplicity, however, we approximate the self energy 
in this investigation in a phenomenological way by the 
single particle energies of the form 



e(fc) = V(fi) + ^/k 2 + M* 2 (n), 



(11) 



where the effective mass M* = AI + S = M+g a a and the 
effective vector field V — g^ui are given by the effective 
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FIG. 2: Diagonal matrix elements (upper part) v(k,k) and 
non-diagonal matrix elements (lower part) v(k, k') with k' = 
0.009 fm -1 as a function of the momentum k in the 1 5o chan- 
nel for Gogny force and the corresponding separable forces. 
The dotted curve represents V| ow k . 



meson fields. They are determined by the solution of 
the nonlinear Klein Gordon equations of the RMF model 
with the parameter set NL3 [55j |. 

First we solve the gap equation Eq. © in the 1 Sq chan- 
nel for nuclear matter at various densities with the Gogny 
force [38J. In order to determine a separable force, we 
first choose for the function p(k) the Gaussian ansatz 



P(k) 



-a A k A 



(12) 



The two parameters G and a are fitted to the function 
A(kp). Comparing with the Gogny force, we obtain two 
sets of parameters G = 738 MeV-fm 3 and a — 0.636 fm 
for the parameter set Dl [| and G = 728 MeV-fm 3 and 
a = 0.644 fm for the set D1S [56$ . In the upper part of 
Fig. [T] we show the pairing gaps in the 1 Sq channel ob- 
tained from the full Gogny forces and the corresponding 
separable approximations (dashed curves). The separa- 
ble forces reproduce the curves for the gap almost per- 
fectly, especially in the case of the parameterization D1S. 
As we see the pairing gap can be very well reproduced at 
low densities by a Gaussian shape and we observe only 
slight deviations for higher densities with kp > 1.0 fm -1 . 
Of course, we could have chosen also some other form of 
p(k) to map the curves of the gap completely. In the 



lower part of Fig. [T] we obtain perfect fits to the results 
of the Gogny force by using a linear combination of three 
Gaussians with different widths and amplitudes for p{k). 



p{k) 



1 3 



(13) 



with the parameters G = 731.25 MeV-fm 3 and a x = 0.69 
fm, <Z2 = 0.16 fm, 0,3 = 0.47 fm for the parameter set 
Dl [H and G = 742.09 MeV-fm 3 and a x = 0.74 fm, an = 
0.26 fm, 03 = 0.36 fm for the parameter set D1S [59 ]. 
In the following calculations, however, we apply only one 
Gaussian form for p{k) that is much simpler and easier 
to handle and that is good enough to describe the pairing 
properties in the finite nuclei. 

In order to investigate the behavior of the separable 
forces in more detail we also plot in Fig. [5] the diagonal 
and non-diagonal matrix elements for each of these cases. 
The matrix elements in the 1 Sq channel of the separable 
forces obtained from Gogny Dl and D1S are very close to 
those of the original Go gny forces. Similar results have 



been found in Refs. [48|, |57J. The Gogny force is an ef- 
fective force in the nuclear medium. In the case of bare 
nucleon-nucleon forces with a strong short range repul- 
sion this would be no longer the case. It is known (58| 
that, by integrating out the high momentum components 
of the bare force, one obtains from most of the realis- 
tic iW-potentials the same interaction a low momenta 
V[ owk . This low momentum NN interaction V la „ k can 
describe the two-nucleon system at low energy very well. 
Although the matrix elements of various realistic bare 
NN interactions are scattered, their low momentum part 
V[ ow k has the same shape. It has been found that the ma- 
trix element of the separable form of AVI 8 is very close 
to that of V 1o „ k obtained from AV18 and very different 
from that of the potential AV18 itself. Therefore we plot 
in Fig. [2] also the matrix elements of V lo „ k . This illus- 
trates the physical content of the separable force as an 
effective interaction for the area of low energies. It is also 
shown in Fig. [2] that the matrix elements of the Gogny 
force D1S have a very similar behavior to that of V[ ow k . 
This indicates that the Gogny force has a clear link to 
the bare NN force, especially for the parameterization 
Dl. 

Now we turn to the solution of the RHB equations in 
finite nuclei. As discussed earlier the separable force has 
simple Gaussian form and it is fitted to reproduce the 
density dependence of the gap at the Fermi surface in 
nuclear matter derived from the Gogny force. This force 
is definitely not identical to the full Gogny force in the 
1 S'o-channel. Therefore, in order to apply these separable 
pairing force in calculations of finite nuclei, one has to 
study to what degree such a separable force can describe 
the pairing properties of these systems. For that purpose 
we use RHB theory in several spherical isotope chains, 
e.g. in Sn- and in Pb-isotopcs. 

First, we transform the force (El from the momentum 
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space to the coordinate space and obtain 

Viv^v'^) = — G (S(R-R') P{r)P{r') |(1 - P CT ) 

(14) 

where R = i(ri + r 2 ) and r = r x — r 2 are center of 
mass and relative coordinates respectively, and P(r) is 
obtained from the Fourier transform of p(k). Using the 
Gaussian ansatz (11211 we find 



(15) 



P(r) 



(47T0 2 ) 3 / 2 



Because of the <5-term in Eq. (|T4|) that insures transla- 
tional invariance this force is not completely separable 
in coordinate space. However, the matrix elements of 
this force can be represented by a sum of a few sepa- 
rable terms in a basis of spherical harmonic oscillator 
functions. 

In order to show this we start from the basis 



\nljrrij) = <p n ij mi (r) = Rni(r, b)[¥i(f) ® xi 



(16) 



where R n i{r,b) = b lR nl (r/b) and [Y t (r) ® x\ ]j mj rep- 
resents the wave function in spin and angels coupled to 
total angular momentum jrrij. The radial wave function 
has the form 
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(17) 



with the radial quantum number n = 0, 1, . . . and the or- 



bital angular momentum I. The quantity b — y/h/(muio) 
is the harmonic oscillator length. In the pairing channel 
we need only the two-particle wave functions coupled to 
angular momentum J — and the projector |(1 — P a ) 
restricts us to the quantum numbers S = L = 0. Recou- 
pling from the LS- to the jj-scheme therefore leads to 
the two-particle wave function 



|12)o = |Vn 1 i 1 j 1 (ri) ! <^„ 2 j 2 j 2 (r 2 ))j=o 



(18) 



= -,R nih {ri,b)R n2l2 (r 2 ,b)\X = 0)\S = 0) 
si 

with j = \J2j + 1 and s = The functions |A = 0) = 
[Y h [h) ®Y h (f 2 )] Q and \S = 0) = [xi ® Xijo are the 
angular and spin wave functions coupled to angular mo- 
mentum A = and spin S = 0. All these wave functions 
are expressed in laboratory coordinates, while the sepa- 
rable pairing interaction in Eq. (| 14[) is expressed in the 
center of mass frame by the center of mass coordinate 
R and the relative coordinates r of a pair. Therefore 
we transform to the center of mass frame using Talmi- 
Moshinsky brackets [sH, [13, HH . We use the definition of 
Baranger [62j 



NLnl 



FIG. 3: Results of RHB calculations in finite nuclei for Sn iso- 
topes (upper panel) and for Pb isotopes (lower panel) : pairing 
energies for the Gogny forces Dl and D1S (solid curves) and 
their corresponding separable forms (dashed curves). 



where 

AffiU = (tfL,nJ,A|n 1 Ii ) n a l 2 ,A) (20) 

are the Talmi-Moshinsky brackets with the selection rule 

2N + L + 2n + l = 2n 1 +h + 2n 2 + l 2 . (21) 

Here we need these bracket only for the case A = 0. We 
therefore can express the wave function (jTSJ) in terms of 
center of mass and relative coordinates 
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(22) 
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The oscillator pa- 



rameters for the center of mass and the relative coordi- 
nates are bjz = 6/v2 and b r — by/2. Finally we find the 
pairing matrix elements of the interaction (|14[) 

V° 2V2 , = (n 1 l 1 j 1 ,n 2 hj2\V\n v l 1/ j v ,n 2 >hj2>)j=o (23) 

as a sum over the quantum numbers N, L, N', L',n, I, n', 
and I' in Eq. ([IT))) . The integration over the center of mass 
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FIG. 4: Pairing energies of Sn- and Pb-isotopes obtained with 
different numbers of separable terms iVo in Eq. (|25[) . 



coordinates R and R' leads to N = N', L = L'. Further 
restrictions occur through the fact that the sum contains 
integrals over the relative coordinates of the form 



Rni{r, b r )Y lm (f)P(r)d 3 r. 



(24) 



They vanish for I ^ and this leads to L = I = 0. 
The quantum numbers n and n' are determined by the 
selection rule (|2"Tj) and we are left with a single sum of 
separable terms 



(25) 



N 



with the single particle matrix elements V{%. For l\ 
h = I, ji = h = j we find 



V& = M%t£A / R n0 (r, b r )P(r)r 2 dr. (26) 
si J 
o 

For a Gaussian ansatz of P(r) in Eq. (fTS")) this integral 
can be evaluated analytically 



N 1 2 1 / 4 (l-a 2 )" j NQn0 V(2n + 1)! 

12 b 3/2 ,,3/4 Q + a 2)«+3/2 g ^l^l 2 n + 1 n\ 

where the parameter a = a/b characterizes the width of 
the function p{r) in units of the oscillator length b and n 
is given by the selection rule (f2"Tj) n = n\ + + I — N. 



Thus we find, that the pairing matrix elements for 
the separable pairing interactions used in the RHB equa- 
tion can be evaluated by the sum of separable terms in 
Eq. (|25[) . In order to study the pairing properties in finite 
nuclei, we solve the RHB equation 
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(28) 



self-consistently for the Dirac Hamiltonian ho and the 
pairing field 



A 12 ^Gj^P, 



N v 12 7 



with the parameters 
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and the pairing tensor k — UV T . The pairing energy in 
the nuclear ground state is given by 



-Epair — G ^ P^P N . 

N 



(31) 



All the following calculations are carried out for the pa- 
rameter set NL3 [55| by expanding the Dirac- Bogoliubov 
spinors in terms of 20 major oscillator shells [20] . 

In Fig. [3] we show the dependence of the pairing energy 
on the neutron number for the chain of the isotopes 100 Sn 
~ 160 Sn and 164 Pb ~ 264 Pb. As we see from the upper 
part of Fig.[3]good agreement is observed for the pairing 
energies calculated with the Gogny pairing force and its 
separable approximation. The largest discrepancy is less 
than 10%. From this comparison we can conclude that 
the separable pairing interaction can describe the paring 
properties of finite nuclei on almost the same footing as 
its corresponding pairing interaction. Therefore, we can 
use the separable pairing interaction instead of its com- 
plicated original form. In fact, one could get even better 
agreement by using the ansatz (|13p with three separa- 
ble terms which produces in momentum space identical 
results as the Gogny force (see Fig.[Tp). 

As we see from the Eq. ([2"5|l the separable pairing in- 
teraction is not fully separable in the spherical harmonic 
oscillator basis. We have a sum over the quantum number 
iV characterizing the major shells of the harmonic oscilla- 
tor in the center of mass coordinate. In practical applica- 
tions it turns out that this sum can be restricted to finite 
values N < Nq — 8. It is therefore enough to determine 
the iVo + 1 matrices at the beginning of the iteration 
and to re-calculate the quantities Pjy m Eq- ([29[) in each 
step of the iteration. As compared to calculations with 
the full Gogny force in the pairing channel this means a 
considerable reduction in memory and computer time. 

In order to study the convergence with the number of 
separable terms iVo we show in Fig. 3] the pairing energies 
in the isotopic chains with Z = 50 (Sn) and Z — 82 (Pb) 
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used for the zero range forces. It has been smoothed by 
a Fermi function 



(MeV) 



-2 2 4 
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FIG. 5: The pairing matrix elements and matrix elements of 
the pairing potential A for the nucleus of 244 Pb. The results 
for various values of No are compared with those for iVo = 
21 which are identical those obtained with the full sum in 
Eq. ((251). 



for various values of No. We find that for the nuclei 
around the line of /3-stability, iVo = 5 is large enough to 
get the full pairing energy; but for the nuclei far from the 
/^-stability line, we need at least N — 8. 

The left panels of Fig. [5] show the size of the individ- 
ual matrix elements (|25l) calculated with a value Nq = 5 
and No = 8 as a function of their exact value (which is 
identical to the value at N = 21). We see that in par- 
ticular the large matrix elements are very concentrated 
along the 45° line. Only for the very small pairing ma- 
trix elements we find small deviations. Similar results are 
obtained for the matrix elements of the pairing field A12 
of the nucleus 244 Pb. Here we find some deviations for 
No = 5. However No = 8 gives satisfactory agreement. 

Finally we have carried out a comparison of the new 
separable pairing force in Eq. (f2"S")) with the full Brink- 
Booker part of the Gogny force in the pairing channel 
(called full Gogny force in the following), and with var- 
ious zero range forces. In the upper panels of Fig. [5] we 
show the size of the corresponding pairing matrix ele- 
ments as a function of the matrix elements of the full 
Gogny force with the parameter set D1S [56| and the 
lower panels present matrix elements of the pairing field 
A12 in the oscillator basis for the finite nucleus 244 Pb as 
a function of those calculated with the full Gogny force 
D1S. In the two panels on the left side we have used the 
separable pairing force with No = 21 which is identical 
to the full sum of Eq. (|2"5)) . The rest of the panels in 
Fig. [6] show calculations with three zero range forces ([I]) 
of various strength parameters. They have been deter- 
mined in calculation for the finite nucleus 244 Pb shown in 
the lower two panels, where a cut-off energy E c has been 



1 



l + exp((£-£; c 



)/D c ) 



(32) 



as in Ref. [29(. The smoothing parameter has been kept 
constant D c = 0.5 MeV in all cases and three different 
values have been chosen for the cut-off energy, E c = 18 
MeV in the lower middle panel and E c = 9, 27 MeV in 
the lower right panel. For each cut-off energy the value of 
the strength parameter Vo has been chosen in such a way 
that the resulting average paring gap in the canonical 
basis 



A Q „ = i£|A 



(33) 



is equal to A a „ = 1.76 MeV, the corresponding average 
gap calculated with the full Gogny force D1S used in the 
abscissa of the three lower curves. We obtain the values 
Vo = 485,326, and 280 MeV-fm 3 for E c = 9, 18, and 27 
MeV respectively. 

We find that although the 5-force can give the same 
average gap as the Gogny force D1S if the size of the 
strength is adjusted properly, the individual matrix el- 
ements of the forces and the matrix elements A of the 
pairing field are very different from each other. Apart 
from many rather small matrix elements the rest of the 
matrix elements of the <5-force cluster around rather con- 
stant values. Therefore the 5-force behaves very much 
like a constant pairing force with a plateau. As seen 
in the two right lower panels of Fig. [5] the value of this 
plateau depends on the cut-off energy. On the other side 
our separable force concentrates along the 45° line espe- 
cially for the large matrix elements. There are only small 
differences observed in the region of small pairing matrix 
elements. For the pairing potential A we obtain a similar 
results in the lower panel of Fig. [5] Here it is clearly seen 
that the separable approximation is very similar to the 
full Gogny force. 

Summarizing, we discuss in this investigation a very 
simple effective pairing interaction in the ^o-channel, 
which is of finite range, translational invariant and sepa- 
rable. This simple force can be easily applied in realistic 
applications of modern relativistic and non-relativistic 
density functional theory, in particular also in compli- 
cated calculations, such as for nuclei with triaxial shapes, 
for nuclei in the rotating frame, for the fission process, 
for QRPA calculations and for all kinds of investigations 
beyond mean field theory using techniques of projection, 
generator coordinates, or particle vibrational coupling. 
Investigations in this direction are in progress. 
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